Physical Structure of the Planetary Nebula NGC 3242 from the 
Hot Bubble to the Nebular Envelope 1 
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ABSTRACT 

One key feature of the interacting stellar winds model of the formation of planetary nebulae 
(PNe) is the presence of shock-heated stellar wind confined in the central cavities of PNe. This 
so-called hot bubble should be detectable in X-rays. Here we present XMM-Newton observations 
of NGC 3242, a multiple-shell PN whose shell morphology is consistent with the interacting stellar 
winds model. Diffuse X-ray emission is detected within its inner shell with a plasma temperature 
^2.35 x 10 6 K and an intrinsic X-ray luminosity ^2x 10 30 ergs s _1 at the adopted distance of 0.55 
kpc. The observed X-ray temperature and luminosity are in agreement with "ad-hoc" predictions 
of models including heat conduction. However, the chemical abundances of the X-ray-emitting 
plasma seem to imply little evaporation of cold material into the hot bubble, whereas the thermal 
pressure of the hot gas is unlikely to drive the nebular expansion as it is lower than that of the 
inner shell rim. These inconsistencies are compounded by the apparent large filling factor of the 
hot gas within the central cavity of NGC 3242. 



Subject headings: planetary nebulae: individual (NGC 3242) 



1. Introduction 

Planetary Nebulae (PNe) consist of stellar ma- 
terial ejected by low- and intermediate-mass stars 
(0.8-1.0 M < Mi <8-10 M ). Towards the 
end of the Asymptotic Giant Branch (AGB), these 
stars experience copious mass loss and eject most 
of their stellar envelope through a slow, dense 
AGB wind. The ejected material is subsequently 
ionized by the central star and becomes a PN. PNe 
eventually disperse into the interstellar medium 
as they expand, whereas the stellar cores, mainly 
composed of carbon and oxygen, will evolve to- 



ward the white dwarf stage. 

Near the time when the hot stellar core is ex- 
posed, the slow AG B wind, with term inal veloc- 
ities 5-30 km s _1 (jEder et al.l Il988l ). is super- 



seded by a fast stellar wind with terminal veloci- 
ties 1,000-4,000 km s" 1 (|Cerruti-Sola fc Perinotto 



1 Based on observations obtained with XMM-Newton, 
an ESA science mission with instruments and contributions 
directly funded by ESA Member States and NASA. 



1985HGuerrero et al.ll2010h . This fast stellar wind 
sweeps up th e slower AGB wind to form a PN 
([Kwokl Il983 ). In this interacting stellar winds 
(ISW) model, the physical structure of a PN would 
be simi lar to that of a wind-blown bubble, as mod- 
eled bv I Weaver et al.l (|1977l ). comprising a central 
cavity filled with shocked fast wind (the so-called 
hot bubble), a dense shell of swept-up AGB wind 
at 10 4 K, and an outer envelope of unperturbed 
expanding AGB wind. In a simplistic model, the 



temperature of the shocked stellar wind inside the 
hot bubble would be 10 7 -10 8 K, but turbulent 



mixing ( e.g.. iMellema fc Frank! 19951) or heat con- 
ducti on (jZhekov fc Perinottolll99a ISteffen et al. 
20081 ) lowers the temperature of the hot gas to 



10 6 -10 7 K and raises its density to produce op- 
timal conditions for the emission of soft X-rays. 
Therefore, X-ray observations of shock-heated hot 
gas in PNe provide us a direct means to exam- 
ine the interaction of the fast stellar wind with 
the AGB wind and to investigate the transfer of 
energy and momentum to the PN envelope. 

ROSAT observations sh owed hints of diffuse X - 
ray emission in a few PNe ( Guerrero et al.l 120001 ). 



However, it was not until the advent of Chandra 
and XMM-Newton, with their unprecedented reso- 
lution and sensitivity, that we were finally able to 
unambiguously detect hot gas in PNe. Chandra 
and XMM-Newton have resolved the diffuse X-ray 

I ■ **-! 

emission in a handful of PNe (e.g., Kastnc r et al 
2000, 120011 IChu et al.|[200lt IGuerrero et al.ll200 ' 



2005) and revealed unexpected, hard X-ray emis- 



sion from the central stars of several PNe that may 
originate from the coronal emission of unseen faint 
binary companions or shocks within the fast stellar 



winds (Guer rero et al.ll2001t iKastner et al.l 12003; 



iMontez et al.112010 ). Observations of diffuse X-ray 
emission in PNe demonstrate that hot gas in ellip- 
tical PNe is confined within the innermost nebular 
shell and that the high pressure of the hot gas may 
indeed drive the nebular expansion as expected in 
bubble models. 

The observed values of Lx and Tx are in agree- 
ment with predic tions of the time-dep endent mod- 
els developed bv lSteffen et al.l (|2008J ) that include 
heat-conduction for PNe with central stars of nor- 
mal, hydrogen-rich surface composition. There 
is, however, an unsolved puzzle: the analyses of 
the chemical composition of the X-ray-emitting 
plasma in PNe suggest that it is mainly composed 
of stellar wind, with little contamination of mate- 
rial from the cold nebular shell (e.g., NGC 2392 



and NGC 6543, IGuerrero et al.1 120051: IChu et al 
20011 ). This discrepancy has been further illus- 



trated by the analysis of the high-spectral reso- 
lution Chandra LETG spectrum of BD+30°3639 
that conclusively confirms the stellar w ind com- 
posit ion of its X-ray-emitting plasma (|Yu et al 
2009). This problem has prompted alternative 



including the action of fast collimated outflows 
and/or sl ow fragments in the ons et of the fast stel- 
lar wind (jSoker fc Kastnerll2003l ) , and the absorp- 
tion of energy from the stellar wind by slowly mov- 
ing ions emb edded in the wind itself, the so-called 
pick-up ions (|Soker et al.ll2010l ). 

NGC 3242 (PNG261.0+32.0), the Ghost of 
Jupiter, is a multiple-shell PN with a bright, 
28"x20" inner ellipsoidal shell and ansae sur- 
rounded by a fainter, 46"x40" moderately el- 
liptical envelope. These two shells are further 
enclosed by arcs and a giant broken halo re - 
vealed by deep images (|Corradi et al.ll2003l 120041 ) . 
The double-shell morphology of the main neb- 
ula of NGC 3242 is highly suggestive of interac- 
tions between the fast ste llar wind of its centra l 
star (^=2,400 km s"\ IPauldrach et al.l l2004h 
and the previous slow AGB wind. The AGB 
wind (the nebular envelope) has been swept by 
the fast stellar wind to form a thin ionized shell 
with a central cavity that can be expected to 
be filled with shocked fast wind. This shock- 
heated gas should emit X-rays, and the diffuse X- 
ray emission from NGC 3242 is likely detectable 
beca use of its proximity (distance = 0.55±0.23 
kpc, iTerzianl fl997l; iMellemal |2004|), and low ex- 
tinction dBalick et al.l Il993l iHenrv et al.l l2000t 



mechanisms for the production of hot gas in PNe, 



Pottasch fc Bernard-Saladl2008l ). 

In this paper, we present XMM-Newton obser- 
vations of NGC 3242 that have detected diffuse X- 
ray emission within its innermost nebular shell. 
The observations are described in §2, the results 
are presented in §3, the physical structure of the 
optical nebula is investigated in §4, and the ef- 
fects of the shocked stellar wind in the nebula are 
discussed in §5. 

2. Observations 

2.1. XMM-Newton X-Ray Observations 

NGC 3242 was observed with the XMM-Newton 
Observatory in Revolution 730 on 2003 December 
4 using the EPIC-MOS1, EPIC-MOS2, and EPIC- 
pn CCD cameras (OBSID = 0200240401). The 
two EPIC-MOS cameras were operated in the Full- 
Frame Mode for a total exposure time of 19.1 ks, 
while the EPIC-pn camera was operated in the Ex- 
tended Full Frame Mode for a total exposure time 
of 15.7 ks. The Medium filter was used for all 
observations. The XMM-Newton products were 



processed using the XMM-Newton Science Analy- 
sis Software (SAS version 10.0.0) and the calibra- 
tion files from the Calibration Access Layer avail- 
able on 2010 September 16. The event files were 
screened to eliminate events due to charged parti- 
cles or associated with periods of high background. 
For the EPIC-MOS observations, only events with 
CCD patterns 0-12 were selected; for the EPIC- 
pn observation, only events with CCD pattern 
(single pixel events) were selected. Time intervals 
of high background, when the count rate in the 
background dominated 10-12 keV energy range is 
>0.3 cnts s" 1 for EPIC-MOS and >1.4 cnts s" 1 
for EPIC-pn, were discarded. The resulting net 
exposure times are 18.7 ks, 18.7 ks, and 13.6 ks 
for the EPIC-MOS1, EPIC-MOS2, and EPIC-pn 
observations, respectively. 

The XMM-Newton EPIC observations detect a 
source of diffuse X-ray emission at the location of 
NGC3242. An inspection of EPIC-pn and EPIC- 
MOS images at different energy ranges reveals 
that this source is soft, with most emission below 
1.0 keV and very little emission at higher ener- 
gies. The EPIC-pn background-subtracted count 
rate in the 0.38-2.0 keV energy range is 31.3±1.6 
cnts ks" 1 for a total of 422±22 counts. The EPIC- 
pn net count rates in the energy ranges 0.38-1.0 
keV and 1.0-2.0 keV are 30.3±1.6 cnts ks" 1 , and 
1.1±0.4 cnts ks" 1 , respectively. The EPIC-MOS 
background-subtracted count rates in the 0.38-2.0 
keV energy band are 4.4±0.5 cnts ks -1 for MOS1 
and 4.8±0.5 cnts ks" 1 for MOS2, with a total of 
83±10 cnts for MOS1 and 90±10 cnts for MOS2. 

2.2. Archival Narrow-band HST Imaging 

In order to examine the spatial correlation be- 
tween the diffuse X-ray emission and the optical 
nebula, narrow-band WFPC2 images of NGC 3242 
in the Ha, He n A4686, [N n] A6583, and [O in] 
A5007 emission lines were retrieved from the HST 
archive (Proposal ID 7501 and 8773, PI: Arsen 
Hajian, and Proposal ID 6117, PI: Bruce Balick). 
The images used in this work are listed in Ta- 
bic Q] with their integration times and filters. In all 
cases, the innermost shell of the nebula was regis- 
tered on the WFPC2-PC1 CCD, while a fraction 
of the outer envelope was missed by the detec- 
tor. These images were calibrated via the pipeline 
procedure and cosmic rays were removed by com- 
bining different exposures obtained with the same 



filter using standard IRAF|j routines. The final 
images have total exposure times of 100 s for Ha, 
160 s for He II, 1260 s for [O III], and 2840 s for 
[N II]. 

2.3. Medium-Dispersion Spectroscopy Ob- 
servations 

Medium-dispersion long-slit spectroscopic ob- 
servations of NGC 3242, obtained using the Gold- 
Cam CCD Spectrograph (GCCAM) on the 2.1m 
telescope at the Kitt Peak National Observatory 
(KPNO) on 1996 December 7, were kindly pro - 
vided to us by Dr. K. Kwitter (JHenrv et aljfeoOCl) . 
The observations, covering the spectral region 
3600-9600 A, were obtained using a 285" x 5" slit 
oriented along the East- West direction and cen- 
tered at 8" south of the central star. The spatial 
scale of the observations is 0'.'78 pix -1 . The grat- 
ing 240 was used with the GG-345 blocking filter 
to cover the blue spectral region 3650 — 6750 A 
at a spectral dispersion 1.49 A pix -1 , whereas the 
red spectral region 5650 — 9600 A was observed at 
a spectral dispersion 1.86 A pix -1 using grating 
58 with the OG-530 blocking filter. Two expo- 
sures of 60 s each were acquired in the blue spec- 
tral range, while exposures of 60 s, 30 s, and two 
of 20 s were acquired in the red spectral range. 
The original spectra were reduced and analyzed 
using standard IRAF routines. For the flux cali- 
bration we used several observations of the spec- 
trophotometric standard stars BD+28°4211 and 
Hiltner600 obtained on the same night. 

3. XMM-Newton Data Analysis 

3.1. Spatial Properties of the X-ray Emis- 
sion from NGC 3242 

To study the spatial distribution of the X-ray 
emission from NGC 3242, we have produced EPIC 
images of NGC 3242 in the energy band 0.25-2.5 
keV by extracting the individual EPIC-pn, EPIC- 
MOS1, and EPIC-MOS2 images, mosaicing them 
together, applying the exposure map correction, 
and smoothing it. The raw image is shown in the 



2 IRAF, the Image Reduction and Analysis Facility, is dis- 
tributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Re- 
search in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 



left panel of Figure Q] while the exposure map cor- 
rected, smoothed version of the image overlaid by 
X-ray contours is shown in the central panel. 

To compare the relative spatial distribution of 
the X-ray-emitting gas and the ionized nebular 
material, we have superposed the X-ray contours 
on the HST WFPC2 [O in] image of NGC 3242 
(Fig. fright). The comparison between the X- 
ray contours and the optical image suggests that 
the X-ray-emitting gas is confined within the in- 
nermost shell of NGC 3242, as the location of 
the lowest intensity level X-ray contours outside 
the 28"x20" inner shell of NGC 3242 is most 
likely caused by the point-spread function (PSF) 
of XMM-Newton that is known to have a half en- 
ergy width (HEW) of 15'/ 2 for EPIC-pn and 13'/0- 
13'.'8 for EPIC-MOS and significantly extended 
wings. 

The X-ray image and contours also suggest that 
the detailed morphology of the X-ray emission is 
asymmetric, with an emission peak northwest of 
the central star. We further compare in Figure [5] 
the surface brightness profile of the X-ray emission 
of NGC 3242 along the major and minor nebular 
axes with the PSF from a point source in the field 
of view. Both profiles of the nebular emission are 
asymmetric and more extended than those of the 
point-source, more clearly for the profile along the 
major axis of the nebula. 

In order to better understand the distribution 
of X-ray-emitting gas in NGC 3242 we have con- 
structed simulated observations for comparison 
with the XMM-Newton images. For these simu- 
lated observations we modeled the emission from 
NGC 3242 as though it arose from a constant- 
density ellipsoidal shell of X-ray-emitting gas inte- 
rior to the innermost nebular shell. We assumed a 
prolate ellipsoid with semi-major and semi-minor 
axes of 11'.'25 x 7'.'5 in the plane of the sky. Three 
models were considered: Model A assumed the 
X-ray-emitting gas was confined to an ellipsoidal 
shell with a fractional width (Ar/r) of 0.1; Model 
B assumed a similar shell with a fractional width 
of 0.2; and Model C assumed the X-ray-emitting 
gas filled the ellipsoid with constant density. 

To simulate yhe XMM-Newton observations, 
we generated random simulated X-rays based on 
the model distribution of the X-ray emitting gas 
with further randomization in the plane of the sky 



consistent with the XMM-Newton PSlH. We also 
added random X-ray events to the simulated im- 
age to mimic the background emission. The total 
number of simulated source and noise counts were 
set to match those of the actual XMM-Newton ob- 
servations. Each model was used to make ten sim- 
ulations to explore the variations caused by the 
small number statistics of these Monte Carlo sim- 
ulations. We then adaptively smoothed the simu- 
lated observations with the same parameters used 
for the actual observations. Figure [3] shows a typ- 
ical realization for each of the three models de- 
scribed. In the realizations of models A and B (the 
shell models) a central deficit in X-ray emission is 
always apparent. In the realizations of Model C 
the X-ray emission peaks on or near the center; 
the example shown in Fig. [3] even exhibited an off- 
set peak matching the actual observation. These 
simulated observations demonstrate that the dif- 
fuse X-ray emission from NGC 3242 is more consis- 
tent with that from a central cavity filled with X- 
ray-emitting gas than from a thin ellipsoidal shell. 
They also suggest that the asymmetric distribu- 
tion of the X-ray emission may be spurius due to 
the low count number. 

3.2. Spectral Properties of the X-ray 
Emission from NGC 3242 

To study the spectral properties of the X- 
ray emission from NGC 3242, we have extracted 
its EPIC-pn, EPIC-MOS1, and EPIC-MOS2 
background-subtracted spectra (Figure [4]). Our 
description of the spectral properties of NGC 3242 
will focus on the EPIC-pn spectrum as the num- 
ber of counts in this spectrum is ^5 times larger 
than in the EPIC-MOS1 and EPIC-MOS2 spec- 
tra. The EPIC-pn spectrum of NGC 3242 is soft, 
with most counts below 1.0 keV. The spectrum 
peaks at 0.5-0.6 keV, and then steadily declines 
towards higher energies. This peak is most likely 
due to the He-like O vn triplet at 0.57 keV. 

We shall note that the shape of the EPIC- 
pn spectrum at energies below 0.55 keV is dif- 
ficult to explain. The unbinned spectrum re- 
veals noticeable oscillations in the count rates with 
small number of counts in energy bins at ^450 
eV and ^500 eV, and large number of counts 



3 Based on XMM-Newton ob servations of the bright, soft 
point source Nova LMC1995 llOrio et al.ll2003f) . 
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Fig. 1.— Raw (left) and smoothed (center) XMM-Newton EPIC images of NGC3242 in the 0.25-2.5 keV 
energy band, and (right) HST WFPC2-PC1 [O III] image of NGC3242. The [O III] image is overplotted 
with the X-ray contours derived from the smoothed EPIC X-ray image. Contours correspond to 10cr, 20er, 
50a, 75cr, IOOct, and 150cr above the background level. 




30 -30 

Radius (arcsec) 



Fig. 2.™ Normalized intensity of the EPIC-MOS X-ray surface brightness profile of NGC 3242 (solid line) 
and a point-source in the field of view of the instrument (dotted line) along the minor (left) and major (right) 
axes of the nebula. 



in energy bins between these two energies and 
at ^400 eV. Since these "spectral features" have 
widths smaller than the EPIC-pn spectral reso- 
lution at this energy range (FWHM~90 eV), it 
is unlikely that these spikes are associated with 
real emission lines. It has been reported that 
optical loading can produce the observed effects, 
but this mechanism can be ruled out because the 
Medium Filter used in our observations prevents 



the optical contamination from point sources as 
bright as m y =6-9 mag. (XMM-SOC-CAL-TN- 
0051), while th e central star of NGC 3242 has 
my=10.3 mag ([van Altena et al.lll995l) . We can 
conclude that these oscillations are most likely 
caused by stochastic effects. To mitigate this issue, 
the EPIC-pn spectrum has been binned to have 
at least 25 counts per channel for further spectral 
analysis. Similarly, the EPIC-MOS spectra have 
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Fig. 3. — Realizations of three Monte Carlo simulations of the XMM-Newton images of NGC 3242 for an 
X-ray-emitting ellipsoidal shell of constant density with a shell thickness 10% (model A) and 20% (model 
C), and a filled shell (model C). Contours correspond to 20(7, 40cr, 80cr, and 120cx above the background 
level. 
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Fig. 4.— (top-panel) EPIC-pn (black), 

EPIC-MOS1 (red), and EPIC-MOS2 (green) 
background-subtracted spectra of the diffuse 
emission of NGC 3242. The best-fit joint model is 
shown as a histogram in the corresponding color. 
The inset shows the nitrogen vs. temperature \ 2 
grid plot of the spectral fit where the black, red, 
and green curves represent the 68%, 90%, and 
99% confidence levels, (bottom-panel) Residuals 
of the best-fit joint model to the EPIC-pn (black), 
EPIC-MOS1 (red), and EPIC-MOS2 (green) 
spectra of NGC 3242 shown in the top panel. 



been binned to have at least 15 counts per channel 
for spectral analyses. 

3.3. Spectral Analysis 

For the spectral analysis, we have adopted 
the nebular chemical abundances (He=0.94Heo, 
C=0.78C Q , N=1.61N , O=0.83O Q , Ne=O.75Ne , 
S=0.20S o , and Ar=0.40Ar o ), and foreground 
hydrogen column density (Ah=5x10 20 cm " 2 ) de- 



rived by iPottasch fc Bernard- Salasl (|2008l ). We 
have then modeled the observed EPIC spectra 
using an absorbed APEC optically thin plasma 
emission mod el and adopting the absorption c ross- 
sections from iMorrison fc McCammonl ( 19831 ). 

This model provides a reasonable fit to the 
EPIC spectra of NGC 3242 's diffuse emission with 
a reduced X 2 of 1.60 (=38.4/24) for fcT=0.190±0.009 
keV (~2.2xl0 6 K), although the best-fit model is 
lower than the observed spectrum at energies <0.5 
keV. The fits are not improved by allowing Ah to 
vary and we find that the best-fit values of kT 
and Ah appear to be correlated (Figure [5]) as 
kT = 0.1975 - 20.0xA H , where kT is given in keV 
and Ah in units of 10 20 cm~ 2 . 

Alternatively, we may allow the chemical abun- 
dances of nitrogen and carbon to vary in the 
fits, as these two elements have spectral lines 
that can contribute in the 0.4-0.5 keV energy 
band. Variations of the carbon abundance do 
not produce noticeable changes in the quality of 




N H [10 22 crrr 2 ] 

Fig. 5. — Temperature vs. column density x 2 grid 
plot of the spectral fit with nebular abundances 
and free column density, where the three curves 
represent the 68%, 90%, and 99% confidence lev- 
els, and the cross marks the best fit for the value 
of 5xl0 20 cm~ 2 adopted for JV H . 

the spectral fit, and thus the carbon abundance 
cannot be well constrained. On the other hand, 
changes to the nitrogen abundances produce a sig- 
nificant improvement of the spectral fit, and the 
best-fit model (Figure ij has kT=0.202t° ^H kcV 
(-2.35 xlO 6 K), N=5.3t?;sN0, and a reduced X 2 
of 1.07 (=24.56/23). The limits to the nitrogen 
abundance set by this fit are better illustrated by 
the nitrogen versus temperature \ 2 grid plot of 
the spectral fit in Figure 01 which shows that the 
range of nitrogen abundances exclude the nebular 
abundance of 1.61 N©. The best-fit value for the 
nitrogen abundance implies a N/O ratio of the X- 
ray-emitting gas to be (N/0)x —1-2, —3.3 times 
greater than the N/O ratio of the ionized nebular 
material, (N/0) nc b —0.36. 

For the best-fit model with enhanced nitro- 
gen abundance, we derive an observed flux of 
(4.2t?'i)xl0 -14 ergs cm" 2 s" 1 and an intrinsic X- 
ray luminosity of (7.3JI°'|) x 10 30 d 2 ergs s _1 in the 
0.4-2.0 keV energy band, where d is the distance 
in kpc. The volume emission measure (EM = 
N 2 V, where N is the electron density and V is 
the emitting volume) of this model is — 5xl0 53 d 2 
cm" 3 . The emitting volume of the inner shell of 
NGC 3242 is -2 x 10 52 e d 3 cm 3 , where e is the vol- 
ume filling factor of the X-ray-emitting gas. Thus, 



the electron density of the X-ray-emitting gas is 
— 5e -1 / 2 d 1 / 2 cm" 3 , and the thermal pressure is 
1.6xl0" 9 e" 1 / 2 d 1 / 2 dyne cm" 2 . At a distance of 
0.55 kpc, the X-ray luminosity of NGC 3242 is 
(2.21q; 2 )x10 30 ergs s" 1 , the electron density of 
the X-ray-emitting gas is 4 6" 1 / 2 cm" 3 , and its 
thermal pressure is 1.2xl0" 9 e" 1 / 2 dyne cm" 2 . 

4. Physical Structure of the Optical Shell 
4.1. Bulk Physical Conditions 

The electron temperature and density of the op- 
tical shell of NGC 324 2 can be derived using stan- 
dard techniques (e.g., iGuerrero et al.lll996[) from 
the long-slit intermediate-dispersion spectroscopy 
presented in §2.3. From these two-dimensional 
spectra of NGC 3242, we have extracted one- 
dimensional spectra of representative regions of 
its inner and outer shells, and measured the fluxes 
of the emission lines using the IRAF task splot. 
Table [2] presents line fluxes normalized to an H/3 
flux of 100. The observed fluxes, F, have been 
dereddened using the IRAF task redcorr to derive 
the intrinsic intensity of the line, /: 



I = Fx 10 



-c Hl3 xf x 



(1) 



where ch0 is the logarithmic H/3 extinction con- 
stant computed by comparing the observed value 
of the Ha to H/3 ratio to the expected theo- 
retical value of 2.87 for c ase B recombination 
(|Oster brock fc Ferlandll2006l ). The observed fluxes 
are subsequently corrected using the values of f\ 
co rresponding to th e inter stellar extinction law 
of ISavage fc Mathisl (|l979h . The Ha/H/3 ratio 
measured in the medium dispersion spectra im- 
plies an extinction coefficient c#,g=0.01 for the in- 
ner shell and 0.05 for the outer shell. We note 
that these values of cnp are smaller than the val- 
ues 0.10-0.15 derived bylBalicket all (|l993l ) and 
Pottasch fc Bernard-Salasl (12008), but consis tent 
with the values derived bv lHenrv et al.1 (|2000i) us- 
ing this same dataset. 

The physical conditions for the inner and 
outer shells of NGC 3242 have been derived us- 
ing temperature-sensitive line ratios of [O in], 
[N n], [O n], and [S III], and density-sensitive 
line ratios of [S n] and [Ar iv]. The ions and 
lines used, their ratios, and the values of electron 
temperature, T , and density, N c , are listed in 
Table [3] The results suggest that the inner shell 



density is ~2,200 cm -3 , in excellent agreement 
with the density derived f r om th e [S n] lines by 
Pottasch fc Bernard-Salasl (J2008I ). and the outer 
shell density is ~370 cm -3 . Among the four tem- 
peratures derived from different lines, we adopt 
that from the [O III] lines, 11,900 K for the inner 
shell and 10,400 K for the outer shell, as the [O III] 
temperature diagnostic is the least affected by the 
extinction law and the lines are the brightest. 
This radial decline in temp erature is consisten t 
with the results presented by iBalick et al.1 ([19931 ) . 
Using these values, we derive thermal pressures of 
3.6xl0 -9 dyne cm -2 and 5xl0 -10 dyne cm -2 for 
the inner and outer shells, respectively. 

4.2. Spatially Resolved Physical Condi- 
tions 

The narrow-band HST images and ratio maps 
of NGC3242 in different emission lines (Figure H]) 
display a simple shell morphology consisting of an 
elliptical inner shell and an attached outer shell. 
This shell morphology is only moderately compli- 
cated by ansae and knots best seen in the [N n] im- 
ages (the so-called Fast Low- Ionizati on Emission 
Regions, FLIERS. liTalick et al.l f l998l) . A spatio 



kinem atic study of NGC 3242 by IBalick et al 



(J19871 ) confirms this simple shell structure: the 
inner shell is a nearly round bubble expanding 
at 25-30 km s _1 , while the outer shell is a co- 
expanding envelope filled with material. Thus, 
the Ha image of NGC 3242 can be used to ob- 
tain density and pressure profiles of the nebula for 
comparison with those of the X-ray-emitting gas 
in the central cavity. 

The density profile of NGC 3242 has b een pre- 
viously determined by ISoker et al.l ([19921 ) using a 
ground-based Ha image. We will use the high- 
resolution HST WFPC2 Ha image to determine 
the electron density, N e , and thermal pressure, 
Pth, of the inner shell and the envelope as a func- 
tion of nebular radius. First, we examined this 
image to select a direction that would provide us 
with a clean surface brightness profile. The cut 
along PA=230° is orthogonal to the inner shell 
rim and covers the full extent of the outer shell. 
Thus, we extracted the surface brightness profiles 
in the Ha, He II A4686, [O in] A5007, and [N n] 
A6583 shown in Figure [7] The surface brightness 
profile of the outer shell is very similar in the Ha, 
[O III], and [N n] lines, but the He n emission is 
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Fig. 7.— Ha, He n A4686, [O III] A5007, and [N II] 
A6583 surface brightness profiles of the inner and 
outer shells of NGC 3242 normalized to the bright 
peak at the rim of the inner shell. The profiles 
are derived by averaging the surface brightness in 
radial bins in a sector along PA = 230° with APA 
= 10°. The inset expands this peak to better il- 
lustrate the spatial differences at this location of 
the emission in the different lines. 



diminished most likely because the inner shell is 
optically thick to He + ionizing photons. As for 
the inner shell, all profiles show the bright shell 
rim, but peak at different radial distances: ~7'.'4 
for He II, ~7'/6 for Ha, and ~7'/8 for [O ill]. The 
emission inside this rim drops steeply, but the fila- 
ments projected in the central cavity produce sec- 
ondary peaks in the Ha, [O III], and [N n] lines, 
but not in the He II line. 

The outer shell has been modeled assuming a 
mean electron density ~370 cm -3 and a temper- 
ature of 10,400 K. Figure Weft shows the fit to 
the Ha surface brightness profile of the outer shell 
of NGC 3242 assuming four different radial depen- 
dences of the density: r~ 2 , r _1 , r -1 / 2 , and r -1 / 3 . 
The model surface brightnesses for the r~ 2 and 
r _1 density profiles decrease outwards too rapidly 
compared to the observed surface brightness pro- 
file, and can thus be excluded. The model surface 
brightness for the r -1 / 2 and r -1 / 3 density profiles 
more closely match the observation. 

Using the r -1 / 2 density profiles in the outer 
shell, we proceed to model the Ha and He II sur- 
face brightness profiles at the rim of the inner 
shell by varying the thickness of the shell and ad- 
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Fig. 6. — HST WFPC2 narrow-band images and ratio maps of the central region of NGC 3242. The images 
are displayed on a squared-root scale, while the ratio maps are displayed on a linear scale. 



justing the ionization fraction of Hc ++ . A shell 
with a constant density 2,200 cm -3 , temperature 
of 11,900 K, shell thickness ~15%, and a sharp 
decay in its inner edge produces a good fit to the 
He II profile. The fit to the rim of the shell in the 
Ha profile is also reasonable, but the model profile 
does not include the contributions from filaments 
to the Ha profile (Figure fright). 

The density profile deduced from the model 
comparisons is shown in Figure ®left. The ra- 
dial profile of the thermal pressure (Figure ® right) 
has been constructed assuming a constant electron 
temperature of 11,900 K for the inner shell, and 
10,400 K for the outer shell, as derived from tem- 
perature sensitive optical line ratios. This figure 
shows that the thermal pressure of the inner shell 
is much higher than that of the outer shell, i.e., 
the inner shell is expanding into the envelope. 

We can also compare the thermal pressure 
of the nebular shell to that of the hot interior. 



Assuming that the hot gas is distributed in a 
constant-density shell interior to but in contact 
with the inner nebular shell, the filling factor can 
be expressed as e = 1 — (ri/r ) 3 , where r and r, 
are the outer and inner radii of the hot gas shell (r 
is also the inner radius of the inner nebular shell) . 
As the hot gas density is ^4e -1 / 2 cm~ 3 (d — 0.55 
kpc), we have plotted the hot gas density as a 
function of its inner shell radius r, in Figure [9]- left 
and the thermal pressure in Figure fright. The 
hot gas shell's inner radii corresponding to e of 
0.15, 0.5, and 0.8 are marked. It is apparent that 
the hot gas pressure exceeds that of the nebular 
shell only for small filling factor values, e <0.15, 
thus requiring the hot gas to be concentrated in a 
thin shell with fractional thickness <0.05. 

5. Discussion 

The properties of the X-ray emissio n from 
NGC 3242, based on a preliminary analysis (JRuiz et al 
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Fig. 8.— (left) Ha surface brightness profile of NGC3242 along PA=230° (thick line) and synthetic Ha 
surface brightness profiles of the outer shell (r >8'/9) for different radial dependences of N c . As for the 
inner shell, a thickness of 15% in radius has been assumed (see below), (right) Ha and He n A4686 surface 
brightness profile of NGC 3242 along PA=230° (thick lines) and synthetic surface brightness profiles of the 
inner shell for different values of the shell thickness: 15%, 30%, and 50%. As for the outer shell, a decay of 
Aetxr -1 / 2 has been assumed (see above). 
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Fig. 9. — (left) Electron density, A^, and (right) thermal pressure, P t h, radial profiles of the central cavity 
and the inner shell and envelope of NGC 3242. The density and thermal pressure of the X-ray-emitting gas 
within the central cavity (radius < 7") are shown as a function of the filling factor, e (red curves). The 
values of these physical conditions for filling factors, e, of 0.15, 0.5, and 0.8 are further labeled on the plot 
at the inner radii of the corresponding hypothetical hot gas shell. 



20061 ) . were compared to ID hydrodynamical mod- 
els of PNe that included the effects o f stellar wind 
evolution and heat conduction by ISteffen et al. 
(|2008[ ). They remarked that the round shape of 



the rim of the inner shell of NGC 3242 makes it 
a suitable PN for the comparison with their ID 
simulations. The revised values of different X- 
ray properties of NGC 3242 presented here need 
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to be discussed in the framework of ISteffen et al 
(J2008l )'s model. This revision mostly affects the 
values of the thermal pressure of the hot bubble 
and inner shell rim, whereas the X-ray luminosity 
(scaled to the distance of 0.55 kpc) and plasma 
temperature of NGC 3242 are basica lly the same 

as those used bv lSteffen et~aD <|2008h . 

We note that the models selected bv lSteffen et al. 



(J2008I ) have a thermal pressure of the hot bubble 
that exceeds that of the rim, while our revised val- 
ues indicate the opposite, i.e., the thermal pressure 
of the rim (P t h=3.6xl0 -9 dyne cm -2 ) is higher 



than that of the bubble (Pth=1.6xl0- 9 e -1 / 2 d 1 ' 2 
dyne cm" 2 ). One possible outcome to this issue 
would be a dis tance for NGC 3242 much far ther 
than 0.55 kpc (JTerzianl [l997i iMellemal |2004 ; at 



a distance of 3 kpc, the thermal pressures of the 
hot bubble and rim would be the same, but this 
would make NGC 3242 too large (r = 0.12 pc) and 
evolved (r — 6,000 yrs) and we consider it to be 
unlikely. 

The other possible solution is to consider a low 
value for the filling factor of the X-ray-emitting 
gas in the bubble; a value e=0.15 would make both 
pressures the same. Such a low filling factor is ex- 
pected to produce a noticeable limb-brightening 
morphology that is not observed (Figures [T] and 
[2]). As shown in Fig. [31 hints of a limb-brightened 
morphology are still apparent for a hot gas shell 
with a fractional thickness of 0.2 (i.e., e ~ 0.5) 
even at the XMM-Newton limited spatial resolu- 
tion (~ 15"). It is arguable, however, that ex- 
tinction has reduced the center-to-limb contrast of 
the X-ray emission, thus shifting the emission in- 
wards, as this is an e ffect expected to be noticeable 
( Steffen et al.ll2008T ) even for the small extinction 



(]Vh=5x10 ' cm 2 ) towards the nebula. 

Alternatively, we must consider the dynami- 
cal effects of the photo-ionization of the nebu- 
lar material, which may have overcome those of 
the currently diminished stellar wind of NGC 324 2 
(JKudritzki et all Il997t iTinkler fc Lamersl I200II . 
In this case, the nebular expansion is driven by 
the thermal pressure increase produced in the rim 
by photo-ionization. This may explain the rela- 
tively large thickness of NGC 3242 inner shell rim, 
as compared to that of NGC 6543. We suggest 
that the rim thickness can be used as an indica- 
tor of the relative thermal pressure of hot gas and 
photo-ionized nebular material; a sharp rim would 



imply that the thermal pressure of the hot gas is 
larger than this of the nebula, and thus that hot 
gas itself is present (and detectable) , while a thick 
inner shell rim would imply a lower pressure for 
the hot gas, if present. 

The low plasma temperature is puzzling in 
view of the enhanced N/O ratio of the X-ray- 
cmitting plasma, about three times higher than 
the N/O ratio of the nebular material. Such 
differences in the chemical abundances seem to 
imply that the X-ray-emitting plasma mostly 
consists of shocked stellar wind, with little con- 
tamination of material from the nebula. The 
low plasma temperature, however, needs large 
amounts of material from the cold nebular shell 
to have been incorporated into the hot bubble. 
This problem, similar to that found i n other PNe 
such as BD+30°3639 (lYu et alJl2009h . NGC 2392 
dGuerrero et al.ll2005h . and NGC 6543 dChu et al. 



20011 ). may require mechanisms for the production 
of X-ray emission in PNe other than mixing and 



heat conduction as summarized by iSoker et al 
(|2010h . 

There is, finally, a noticeable asymmetry of 
the spatial distribution of the X-ray emission 
from NGC 3242, with the northwestern half of 
the hot bubble being brighter than its southeast- 
ern half (Figures Q] and [2|). Even if we assume 
that the inner shell of NGC 3242 were a tilted el- 
lipsoid with its Northwestern tip moving towards 
us (and its Southeastern tip receding from us), it 
seems unlikely that an uneven spatial distribution 
of intervening material could be responsible of 
the observed asymmetry as it ma y be the case 
of th e dusty PNe B D+30°3639 (JKastner et al. 



2000h and NGC 7027 dKastner et all l200llh be- 
cause only small ext inction variations were ob- 
served in NGC 3242 (JBalick et al.lll993h . On the 
other hand, our simulations in Section 3.1 suggest 
that the asymmetric distribution of X-ray emis- 
sion could result from small number statistics. 

6. Summary 

We have obtained XMM-Newton X-ray obser- 
vations of NGC 3242, a multiple-shell PN consist- 
ing of an inner shell with a bright, round rim, and 
an outer envelope. The observations have detected 
diffuse, soft X-ray emission confined within the 
innermost shell of NGC 3242. The relatively low 
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temperature of the hot gas, Tx=2.35xl0 6 K, com- 
pared to the expected adiabatic post-shock tem- 
perature _ofja 1 _stenaj^_wind .with a velocity of 2,400 
km s _1 (jPauldrach et al.ll2004f) suggests that heat 
conduction has taken place. Indeed, models in- 
cluding heat conduction provide a reasonable de- 
scription of th e X-ray temperatur e and luminosity 



of NGC3242 (jSteffen et al.ll2008J ). However, the 



chemical abundances of the X-ray-emitting plasma 
are closer to the stellar values, suggesting little 
evaporation of cold nebular material into the hot 
bubble and thus contradicting the expectation of 
heat conduction. 

We have compared the physical properties (N e , 
T e , Pth) of the gas in the hot bubble with those of 
the gas at the inner shell rim. The inner shell can 
be described as a thin shell with a constant density 
of 2,200 cm -3 and a thickness 15% its radius, while 
the envelope is best described by a shell whose 
density declines ex r -1 / 2 to r -1 / 3 . The gas in the 
hot bubble has lower thermal pressure than the 
gas in the shell rim, unless the X-ray-emitting gas 
is mostly confined within a thin shell, which is 
not supported by the X-ray morphology observed. 
Comparisons with simulations favor a large filling 
factor for the X-ray-emitting gas. 

We also note the asymmetric distribution of the 
X-ray emission within the inner shell of NGC 3242 
and found it unlikely to be a result of nonuniform 
absorption of the X-ray emission within the opti- 
cal shell. X-ray observations at higher spatial res- 
olution and signal-to- noise ratio (Montez et al., in 
preparation) are needed to provide a sharper view 
of the distribution of the X-ray-emitting plasma 
within NGC 3242. 
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Table 1 
HST WFPC2 Observations of NGC 3242 



Emission Lines Number of images i C xp Program ID 

M 

Ha 1 

Ho II 1 

[N n] 3 

2 

4 
[O in] 1 

5 



Table 2 
Line Strengths for NGC 3242 



100 


6117 




160 
400 


6117 
6117 




300 


7501, 


8773 


260 


7501, 


8773 


260 


7501 




200 


7501, 


8773 



Inner Shell Outer Shell 



Lino ID Wavelength 

(A) 



[O ii] 3726 0.26 10.8 11.7 3.7 

[On] 



[Sn] 

[Sn] 
H(5 
H7 
[O in] 
[Ar iv] 
[Ar iv] 
H/3 
[O in] 
[O in] 
[CI in] 
[Nil] 
[S in] 
[Nil] 
Ha 
[Nil] 
[S 11] 
[S 11] 
[On] 
[On] 
[S in] 
[S in] 
F(H/3) (erg 



3729 


0.26 


9.7 


10.5 


0.3 


0.3 


4068 


0.20 


0.5 


0.5 






4076 


0.20 


0.5 


0.5 






4101 


0.18 


26.4 


25.6 


27.2 


25.5 


4340 


0.135 


45.5 


47.0 


46.2 


48.2 


4363 


0.13 


10.9 


11.3 


10.1 


10.5 


4711 


0.04 


4.8 


4.9 


4.0 


4.0 


4740 


0.04 


4.0 


4.0 


2.9 


2.9 


4861 


0.00 


100.0 


100.0 


100.0 


100.0 


4959 


-0.02 


339.1 


336.8 


362.2 


359.7 


5007 


-0.03 


1018.1 


1007.8 


1523.3 


1507.9 


5517 


-0.15 


0.2 


0.2 


0.3 


0.3 


5755 


-0.21 


0.07 


0.07 






6312 


-0.30 


0.7 


0.6 


0.4 


0.4 


6548 


-0.34 


0.7 


0.6 






6563 


-0.34 


289 


287 


296 


286 


6583 


-0.34 


2.1 


1.9 


14.6 


13.4 


6716 


-0.36 


0.3 


0.3 


1.1 


1.0 


6731 


-0.36 


0.4 


0.4 


1.0 


0.9 


7320 


-0.43 


0.5 


0.4 


0.05 


0.04 


7330 


-0.43 


0.5 


0.4 


0.05 


0.04 


9069 


-0.64 


3.7 


3.1 


22.6 


19.2 


9532 


-0.65 


18.0 
3.1X10" 11 


15.2 


18.1 
2.2xl0" 12 


15.2 



Table 3 
Physical Conditions in NGC 3242 



Observed Ratios Value 



Physical Parameter Ion Line Ratios Inner Shell Outer Shell Inner Shell Outer Shell 

T e [Nil] (6548+6583/5755 41.6 ... 14,700 K 

T„ [S in] (9069+9532)/6312 32.0 106.6 14,040 K 8,070 K 

T e [O in] (4959+5007)/4363 124.9 189.0 11,880 K 10,400 K 

T e [O n] (3726+3729)/(7320+7330) 19JJ 43_3 10,100 K 9,860 K 

"lVe [Ar IV] 4711/4740 l7l8 l737 2,250 cm" 3 400 cm" 3 

JV e [S n] 6716/6731 0.71 1.13 2,200 cm" 3 340 cm" 3 



14 



